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^ : Abstract 

■ Dissipative dynamics of a quantum bit driven by a strong resonant field and 
. interacting with a heat bath is investigated. We derive generahzed Bloch 

equations and find modifications of the qubit's damping rates caused by Rabi 
\ oscillations. Nonequilibrium decoherence and relaxation of a phase qubit in- 

\ ductively coupled to a LC-circuit is considered as an illustration of the general 

■ results. It is argued that recent experimental results give a clear evidence of 
Q^ ■ effective suppression of decoherence in a strongly driven flux qubit. 
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I. INTRODUCTION 



Many years ago I. Rabi has shown [1] that a population difference of a two-level system 
(or a qubit in a modern language) subjected to the action of a resonant electromagnetic field 
oscillates in time with a frequency proportional to the field strength. Since then the Rabi 
oscillations and corresponding Rabi splitting of energetic levels have been observed in atomic 
ensemble coupled to a single- mode cavity [2], in a single excitonic quantum dot [3-6] as well 
as in a Cooper-pair box [7] and in a Josephson junction [8]. A successful demonstration 
of Rabi oscillations in a Coopcr-pair box combined with Josephson junctions [9] along with 
completed and projected experiments on a coherent manipulation of flux qubits [10-13] offer 
strong possibilities for a practical realization of a superconducting quantum bit [14]. This 
objective can be accomplished by suppressing noise and decoherence caused by a dissipative 
environment (a heat bath). 

Corresponding decoherence and relaxation times, T2 and Ti, are usually expressed in 
terms of the spectral density of the heat bath fluctuations, J(a;), which represents yet an 
imaginary part of a heat bath susceptibility xi^)^ — x"i^) [15-17] 



and T is temperature of the environment. Henceforward the Boltzmann and Planck con- 
stants are supposed to be unity, ^ = 1,^^ = 1, unless otherwise specified. A combination 
S{lj) = J{uj) coth{uj/2T) stands here for the spectrum of heat bath fluctuations. Spectra 
of additional noise sources should be properly added to S{u!) in Eq.(l), in particular, a 
spectrum of 1// (flicker) noise, Sf{u!) ~ l/li^l, which is inherent to every electrical circuit. 
These circuits are necessary to control and read out the state of the superconducting qubit. 
It follows from Eq.(l) that in the presence of energy bias (e 7^ 0) 1//- noise gives an infinite 
contribution to the transversal rate that leads to very fast dephasing of the qubit. Ob- 
viously a formal application of Eqs.(l) to the case oil/ f noise can not be justified because 
the formulas (1) are valid only for the weak interaction between the qubit and the heat 
bath. An additional point to emphasize that the above-mentioned expressions present the 
rates of qubit 's relaxation to the thermodynamically-equilibrium state in the absence of any 
time-dependent external force acting on the qubit. With the driving force the relaxation 
and decoherence rates are modified as it is demonstrated by F. Bloch [18] for a spin 1/2 
in a rotating magnetic field. This phenomenon opens the road for the suppression of deco- 
herence by the external time-dependent field which is especially effective when the specific 
time scales describing the internal qubit 's dynamics are less than the correlation time Tc 
of the heat bath. The similar mechanism underlies a concept of a quantum "bang-bang" 
control, proposed by Viola et al [19,20]. It was shown that a sequence of microwave pulses 
can completely suppress an environment-induced decoherence of the qubit if the repetition 
frequency of the pulses exceeds the inverse correlation time of the heat bath. The gen- 
eral concept is exemplified by an exactly solvable model of a two-level system with energy 
splitting fluctuating due to interaction with a heat bath. 




(1) 



Here A and e are the tunnel splitting and the energy bias of the qubit, t^o = VA^ + e^. 



2 



Effects of the strong driving force on decoherence and relaxation of the dissipative two- 
level system were under extensive discussion for many years [16,21-23]. Bloch-Redficld type 
equations with an arbitrary control field have been derived [24] in the process. However, 
in addition to the general theory there is an urgent need now in the transparent analytical 
formulas for the decoherence and relaxation rates of the resonant-driven quantum bit. 

In the present paper wc analyze efi^ects of Rabi oscillations on decoherence and relaxation 
of a two-level system (a qubit) with the goal to find field-induced modifications of the 
formulas (1). To do that, we resort to the theory of open quantum systems [25-27] as well 
as to the rotating wave approximation. This approximation is quite justified here, because 
for the most experimental situations the frequency of the Rabi oscillations induced by the 
resonant driving force is much less than the energy splitting of the qubit. 

The paper is organized as follows. In Sec. II we derive non-Markovian equations for the 
qubit's variables in the rotating frame and bring them to the form of the generalized Bloch 
equations with kinetic coefficients depending on the Rabi frequency. Dissipative dynamics 
of the qubit is described in Sec. III. We calculate the nonequilibrium damping rates and 
illustrate the results using the fiux qubit coupled to a LC-circuit (a tank) as an example. A 
theoretical disccussion of recent experimental results [12] for dephasing and relaxation times 
of strongly driven fiux qubit is presented in Sec. IV. 



II. GENERALIZED BLOCH EQUATIONS 

The Hamiltonian of the qubit interacting with a resonant field F{t) and coupled to a 
dissipative environment with a variable Q{t) has the form 

H ^^a^ + ^a^-a^^Fcosuot-^a^Q, (2) 

where A is a tunneling rate of the qubit, e is a bias, and Fcosa;ot is a resonant driving 
force with the frequency exactly equal to the energy splitting of the qubit: uq = VA^ + e'^. 
Here we have some distinctions from the Hamiltonian describing a dissipative evolution of 
a spin in a rotating magnetic field [18,23,26]. First of all, in the present formulation of the 
problem the alternating external magnetic field is aligned with one axis z. Because of this 
equations describing a free evolution of the qubit (without its coupling to the heat bath) do 
not have exact solutions and we have to resort to the rotating wave approximation. Besides 
that we have an energy bias e in Eq.(2), and the heat bath variables Q are coupled only to 
z-direction of the qubit's Pauli matrix. For the superconducting fiux qubit [28,29,11] the 
parameters e, F and Q can be considered as static, oscillating and fluctuating parts of the 
flux produced by the external circuit. 

It is convenient to introduce the following operators describing the qubit in the rotating 
frame of reference: 

A = — H cr^, 

UJo Uo 

Y — aycoscuot + — az cr^ smouot, 

\ljJ0 LOo J 

Z= ( — — ^CTa;) COS a;ot — CTy sin a;ot. (3) 
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These operators explicitly depend on time, but preserve all commutation rules inherent to 
Pauli matrices, for example, [X, y]_ = 2iZ, and so on. In the rotating wave approximation 
(RWA) the effective Rabi frequency of the qubit Qr = {A/ujq)F is much less than the 
frequency of the driving force ujo,Qji <^ Uq, so that a parameter a,a — flji/iujo, is much 
less than one, a <S 1. Then, the Hamiltonian (2) of the system has the form: 

H = ei±^x -^Z-U- {Z cos coot + Y sin Cot) + -x\q. (4) 
2 2 2 t cuq coq J 

Here we take into account first-order corrections to the RWA described by the parameter 
a and resulting in the frequency shift of the two-state system. In the absence of the driv- 
ing force and dissipation the Hamiltonian of the qubit is proportional to the operator X. 
Then, the time evolution of diagonal elements of the qubit 's density matrix in the energy 
representation is determined by X, whereas the operators Y and Z determine the evolution 
of the off- diagonal elements of the density matrix. With the driving force the situation is 
essentially complicated. Nevertheless, by convention the decay time of the operator X will 
be subsequently referred to as a relaxation time, and the decay time of the operator Z will 
be referred to dccohcrence time [30]. 

We resort to the rotating wave approximation and drop fast oscillating terms. Then the 
time evolution of the Heisenberg operators X, Y, Z is governed by the following equations: 

X — QrY + — (Y cos ujot — Z sin uot) Q, 

LUq 

Y = -QrX - aQRZ - (—Xcosuot - —z) Q, 



Z = ailnY + (—X sin tuot - —y) Q. (5) 

In the absence of the qubit-bath interaction the variables X, Y oscillate in time, 

X{t) = X{ti) cos Unit - ti) + Y{ti)smnR{t-ti) - aZ{ti)[l - cosilnit - ti)], 
Y{t) = F(ti)cosfiij(t-ti) - X{ti)smnR{t-ti) - aZ(ti) sin - ti), 

Z{t) = Z{ti) + aY{ti) sin QR{t-ti)- aX{ti)[l - cos QR{t - h)] . (6) 

Because of this, the parameter D,r, 

represents the Rabi frequency. The dissipative terms in the equations (5) contain fast 
oscillating components which should be subsequently omitted. In the case of the Gaussian 

fluctuations of the unperturbed heat bath variables or for weak coupling between the 
qubit and the heat bath the total operator Q{t) is a linear functional of the qubit's variables 
[25,31]: 

Q{t) = Q(°)(i) + 1 f dtiip{tM) [— (Z{ti) cosuoh + Y{ti) smuoh) + —X(ti)\ . (8) 
Here 
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^(t,h)^{i[Q^'\t),Q^'\h)U)e(t-h) (9) 

is a linear response function of the heat bath (a retarded Green function). By averaging 
the equations over the initial state of the heat bath with the temperature T we obtain the 
system of the non-Markovian equations: 

A r 

{X) = nR{Y) + — 2 / dti{M{t,ti){i[Y{t) cos (jJot - Z{t)smuJot, 

AZ(ti) cosa;oti + AY(ti) sin(x;oti + eX(ti)]-) + 
(p{t,ti){{l/2)[Y{t) cosujot - Z{t) sinujot, 
AZ{ti) cosa;oii + /SY{ti) sinu;o^i + eX(ii)]+)}, (10) 

(Y) = -nniX) - aQniZ) - ^/ dh{M{t,h){i[AX{t)cosu;ot - eZ{t), 

AZ{ti) cosuJoh + Ar(ti) sincjoti + eX(ti)]_) + 
V?(t,ti)((l/2)[AX(t) cosuJot - eZ{t), 
AZ{ti) cosu;oti + AF(ti) sinc^o^i + eX(ti)]+)}, (11) 

(Z) ^anR{Y) + ^ [ dti{M{t,ti){i[AX{t)smu}Qt-eY{t), 

AZ{ti) cos(x;oti + AY(ti) sinu;oti + eX(ti)]-) + 
ip{t,ti) ((1/2) [AX{t) sin a;oi - eY{t) , 
AZ{ti) cosa;oii + AY{ti) sina;o^i + eX(ti)]+)}. (12) 

Here M(t, ti) = M{t, ti)d{t - ti), and 

M{t,t^) = {{l/2)[Q^'\t),Q^'\t^)U) (13) 

is the symmetrized correlation function of the heat bath fluctuations, h — l,kB — i- Hence- 
forward we use notations [...]_ and [■■■]+ for commutators and anticommutators, respectively, 
and 9{t) is the Heaviside step function. It should be noted that the spectral density of heat 
bath fluctuations, 

S{u) = J dre'^'^Mir), (14) 

is related to the imaginary part x"{'^) = It^xiy^) = Jil^) of the susceptibility 

X(a;) = / dre-Xr) (15) 

according to the fluctuation-dissipation theorem 

5((^)=x"(a;)coth(^). (16) 

To derive the system (10), (11), (12) we resort to the general theory of open quantum systems 
developed in Refs. [25-27,31]. These equations are valid exactly in the case of Gaussian 
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fluctuations of tlic unperturbed variables of the heat bath {Q^^\t)} or approximately for 
the weak coupling between the qubit and the heat bath. In the present paper we restrict 
ourselves to the Bloch-Redfield approximation when the qubit-bath interaction is not strong, 
so that relaxation of the qubit during the correlation time Tc of the heat bath is negligible, 
whereas the period of Rabi oscillations, 2t^/Q.r, can be of order of Tc : ^rTc ~ 1. Generally 
speaking, the perturbation theory is not quite applicable to the case of 1/f spectrum of 
heat bath noise as well as to the case of resonant coupling of the qubit to the high quality 
LC circuit having a peaked spectral density. But even for these cases the perturbative 
approach allows us to find conditions when the contribution of the environment to the 
dissipative dynamics of the qubit will be moderate, thus justifying using the Bloch-Redfield 
approximation. 

The non-Markovian equations (10)-(12) contain the commutators of the qubit's operators 
taken at different moments of time t and ti separated by the correlation time of the heat 
bath Tc- To calculate the commutators we express the operators X(ti),F(ti), and Z{ti) in 
terms of the operators X{t), Y{t), and Z{t) using the free evolution of the "dressed" qubit 
(6). It should be noted that reducing the operators and Z{t) to the operators 

results in the same results for the qubit's damping rates T,j.,Vy,Vz- The 
operators of the qubit X, Y. Z taken at the same moment of time t obey the well-known 
commutation rules for the Pauli matrices, i[X(t), = —2Z{t), and so on. Then, the 

commutators, hke i[X{t),Y{ti)]-, oscillate with the Rabi frequency as functions of the 
time difference t — t — ti : 



The same is true for the anticommutators: {l/2)[X(t),Y(ti)]+ = sinflRT. These 
(anti)commutators should be substituted into Eqs.(10)-(12). Integrating over r gives us 
the Fourier transforms of the correlation function (14) and the response function (15) of the 
dissipative environment taken at the frequencies c^oj^o i and Qr. After averaging the 
collision integrals in Eqs.(10)-(12) over the fast oscillations with the frequencies ujo,2ujo,... 
{cosLUotcosLUoti ~ (1/2) COS a;oT, ...) and omitting unimportant terms we obtain the gener- 
alized Bloch equations for the qubit's variables averaged over the initial state of the heat 
bath: 



Here we drop the brackets (..) symbolizing averaging over the heat bath fluctuations. It 
should be mention also that in order to find the physical values we have to average the 
operators X, Y, Z over the initial state of the qubit. The relaxation coefficients in Eqs.(17) 
are determined by the spectrum of heat bath fluctuations S{u!) (14) taken at the frequencies 



t[X{t),Yit^)] 



2Z{t) cos Qrt + 2aY{t) sin ^rt. 



X + T^X = nRY-u,, 

Y + TyY = -QrX - OQrZ + Uy, 

Z + T^Z = uVLrY + 



(17) 



^0, ^0 i ^R '■ 



A2 r 



(1 + 2a)S{uJo + Qr) + (1 - 2a)S{uJo - Qr) 




A2 r 



S{(jJq) + a 



(18) 



r 



y 
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(1 + a)S{u;o + Qr) + (1 - a)S{u;o - Q r) , e 



+ 



R) 



(19) 



The parameters u^., Uy, defining the steady-state values of the qubit's operators are pro- 
portional to the imaginary and real parts of the heat bath susceptibility x(c<;) (15): 

A2 r (1 + 2a)x^^(a;o + ^r) + (1 - 1oi)yl'{u:^ - ^r) ' 
^^ = 4^ ^ ^""'^^ 2 



[I - q;)x'(^o - ^r) - (1 + «)x'(c^o + ^r) 



+ ax'(^o) 



«^ [x'(o) - xHS^r)] , 



[1 + a)x"(a;o + ^^i?) - (1 - a)x"(c^o - ^^s) 



- ax'{uJQ) 



+ 



X"(f^i?). (20) 



The above-mentioned expressions for the kinetic coefficients of the qubit have been derived 
with taking into account first corrections to the rotating wave approximation. 

III. NONEQUILIBRIUM DISSIPATIVE DYNAMICS OF THE QUBIT 



For the steady-state values of the qubit's operators Xq, Fq, in the rotating frame of 
reference we find the following expressions 



VtRVy 



o2 I r r 

1 Lj^ T J- a; J- y 

^Rl^x + i x^y 

02 _|_ f f 
-r L x'- y 



(21) 



A^[x^^(a;o + Qr) - x"{u;o - ^r)] + 46^X^X^fl) 
A''[S{u;o + nR)+S{u;o-nR)]+^e^S{nR) ' 



(22) 



The small corrections to the rotating-wave approximation have been omitted here. In the 
equilibrium case, when Qr = 0,10 = 0, Zq = 0, we obtain the well-known expression for 
the difference of population Xq,Xq = —Vx/'^x = — tanh(co'o/2T). At the same time, in 
the strong nonequilibrium case with VIr ^ T the steady-state values of X and Y go to 
zero: Xq = 0, Fq = 0, so that the populations of the excited and ground states of the 
qubit are equalized under the action of the resonant field; in so doing the qubit is polarized 
along Z-axis. The dissipativc evolution of the operators X, Y, Z from their initial values 
X(0), y(0), Z(0) is described by the equations 



X{t) = X{0)e''/^' (cosQct 



-t/Ti 



COS Qct 



sin Qct] + 



2Q 



c 



[r(o)-ro] e-*/^^^sinnc^, 



Y{t) = F(0)e-*/^^ (cosl^c^ 



sin flct ] + 



2Qc 
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Xa 



(cos D,ct — sin D,ct\ + 

[X(0)-Xo]e-*/^^^sinQci, 
Z{t) = Z{0)e-^^' + Zo(l- e"^^*) . (23) 

Here 



^^c = y^^|-|(r.-r,)2 (24) 

is the effective Rabi frequency sliiftcd because of relaxation processes, and Ti^, 
r. + r, A2 (1 + 3a)S{u;o + Qn) + (1 - Sa)S{u;o - iln) ^ 



(25) 



is tfie longitudinal damping rate (with non-RWA amendments) which determines the relax- 
ation of the qubit energy. It is worth to mention that the fast oscillating terms omitted 
in the process of our calculation, which is based on the rotating wave approximation, re- 
sult in the corrections to the relaxation rate T^^ (25) of order of the small parameter 
3a = {3nR/Aujo) < 1. 

Under the strong driving conditions, D,r ^ F, the low-frequency noise contributes to the 
relaxation rate Tf^ as well. In the absence of the driving force and the Rabi oscillations, 
when F = 0, Qr = 0, the effective Rabi frequency (24) becomes imaginary, Qc = (V2)(rx ~ 
Ty), and damping of the qubit's variables X, Y, Z is determined by the corresponding rates 
Tx,Ty,Tz, respectively. In this equilibrium case the formulas (1) come into play; in so doing 
the parameter is related to the longitudinal relaxation rate, — Tf ^, whereas Fy and 

describes the dccohcrcnce of the qubit caused by its interaction with the heat bath: 



Tj/ = = T2 ■ Interestingly that with the driving force the dissipative evolution of the 
operators X and Z is also determined by the different time scales Ti (25) and (19). 

Let's assume that in the initial moment of time, t = 0, the qubit was in the ground 
state of the Hamiltonian ijJqX/2. Averaging over this state gives: (-'^(O)) = —1, (^(0)) = 
0, {Z{{))) = 0. Then, with the strong nonequilibrium conditions, ^ F, the probability to 
find the qubit in the excited state, Pe — {1 + {X)) /2, demonstrates damped Rabi oscillations 
with the relaxation time Ti (25) 

PE{t) = ^ (1 - e-*/^^ cos ^Int) . (26) 

At large times, t ^ Ti, the excitation probability Pe tends to 1/2 that points to the 
equalization of level populations. 

The distribution of the qubit between left and right wells can be described by the 
difference of probabilities to find the system in the left well, Pl, and in the right well, 
Pr '■ Pl — Pr — (cz), with the total probability Pl + Pr — 1- Then, for the qubit, ini- 
tially prepared in the ground state, this difference also oscillates in time, but with several 
frequencies: cuq, <^o ± ^r, and flu : 
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{Pl - PrW) = {a,) = — e-*/^i sin sin a;oi + 

—Zq fl - e"^^*) cosuQt - — e-*/^i cosQiji. (27) 

In the initial moment of time we have: -Pl(O) — -Pr(O) = — e/u^o < 0, that means that the 
qubit was shifted towards the right well because of the bias e. Without the bias the differ- 
ence between the probabilities to find the qubit in the left and right wells shows quantum 
beatings with the frequency uj^ corresponding to the energy splitting, the envelope amplitude 
therewith oscillates with the Rabi frequency Qi?: 

Pl-Pr = e~*/^i sin VLrI sin UQt + Zo{l- e~^^*) cos ujot. (28) 

It is of interest that the Rabi oscillations of the envelope are damped out for the time scale 
Ti (25) of energy relaxation, whereas the steady-state polarization of the qubit along the 
direction of the external oscillating field is estabhshed for the different time scale (19). 
In the presence of the bias the nonequilibrium damping rates T^^ and are related as 

For the three-junctions qubit [10,28,29] a persistent current in the qubit loop is proportional 
to cr^— matrix: / = /gCr^, so that the function (P^ — Piijit) = (az) describes flipping of the 
current and the flux in the qubit's loop. 

It follows from Eqs.(19),(25) the Rabi oscillations have a profound impact on decoherence 
caused by 1// (flicker) noise. Contrary to the equilibrium case, in the presence of Rabi 
oscillations the dominant contribution of 1// noise into damping rates of the biased quantum 
bit is finite: (F^)/ = 2{Ti^)f = [e^ /2uJq){v /VLr)^ where Sf{uj) = v/\uj\ is a spectrum of 1/ f 
noise with a frequency-independent constant v. No wonder that we do not have any divergent 
contributions now, because the low-frequency flicker noise is averaged over the rapid Rabi 
oscillations of the qubit. Our pcrturbative approach is valid if the flicker noise constant v is 
smal enough, so that the contribution of 1// noise to the damping rate, (F^)/, is much less 
than the Rabi frequency Qr: {e^ /2uJq){v/Qr) -C 1. 

Besides the clear physical meaning the signiflcance of the expressions (19), (25), (26)- (29) 
lies in the fact that these formulas give us a simple description of the dissipative evolution of 
the various Rabi-oscillating quantum bits. It is important, in particular, for the experiments 
with the phase qubits [10,28,29,13] coupled to a LC-circuit (tank). The resonant tank circuit 
serves as a control element but also as a source of noise with a frequency-dependent spectral 
function S{u) (16). A persistent current in the qubit loop with an inductance Lq can take two 
values ±Iq. The circuit has an inductance Lt, a capacitance Ct and is characterized by the 
resonant frequency cut — 1/\^LtCt and the quality factor Qt = '^tIi- Here 7 = 1/{RtCt) 
is a linewidth of the circuit having a resistance Rt- The qubit and the tank are coupled 
by a mutual inductance = k^LxLq with k'^Qx — 0.1 — 1, so that an interaction term 
has the form: Vint = —McIpIrO'z where It is a current in the tank. It means that the 
heat bath operator Q{t) involved in the Hamiltonian (2) is proportional to the current 
It '■ Q = 2McIqlT- After calculating the susceptibility and the spectrum of the LC-circuit 
we flnd the function x"(a;) (15): 
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x'u = tf^?v; (,._,g^^.,. . (30) 

The resonant character of the heat bath spectrum S{lj), S{lj) = x"{^) coth(ti;/2T), results in 
the dependence of the relaxation and decoherence rates of the qubit, T^^ (25) and (19), on 
the Rabi frequency (7). This dependence is especially pronounced at low temperatures, 
T < ^.R. In particular, in the case when the tunneling splitting cuq = \/ + is much more 
than the resonant frequency of the tank, cuq ^ ujt, the energy distribution Pe (26) and the 
distribution of the qubit between wells Pl — Pr (27) relaxes with the rate 



T,- = kW,^^-^ [4coth (^) + coth (^^^] + coth (^^^) 



+ 



e2 



^'^^^^^^'^(Q|-J^^\7^Q|'°'^ (^) • ^^'^ 

It is known [10,28,29] that for the persistent current qubit the non-zero bias e is necessary for 
the measurements. At the same time the non-zero bias means the significant contribution 
of the last term in Eq.(31) into the relaxation rate Tf"^. This contribution can be especially 
dangerous if the Rabi frequency Qr approaches the resonant frequency of the tank uj^ : Qr ~ 
ujt- To estimate the relaxation time Ti we consider the phase qubit with the inductance Lg = 
lOpH, and the current Iq = 400nyl. Then, the flux created by the qubit's loop is evaluated 
as: Lgig = 2 X 10~^$o, where $o = h/2e is the flux quantum, and Lgl^ = 1.6 x 10~^^J. In 
the case of the exact resonance between the Rabi frequency and the frequency of the tank, 
Qr — ujti the major contribution to the damping rate Tf ^ is given by the simple formula 



rr' = A=<3T^i#^coth(^). (32) 



ul n \2T 

With k'^Qx = 0.1, e/A = 0.1 we find that the Rabi oscillations of the qubit are terminated 
for the time Ti = 60 nsec at low temperature T <^ VLr. The steady-state polarization of 
the qubit is established for the shorter time ~ (l/2)Ti ~ 30 nsec. Under the above- 
mentioned conditions it is possible to observe the Rabi oscillations of the qubit resonantly 
coupled to the tank {VLr = ojt) when, for example, the qubit has the tunneling frequency 
A = 10 GHz, and the Rabi frequency VLr is about 1 GHz. The qubit undergoes many Rabi 
oscillations during the relaxation time, ^rTi ~ 60 ^ 1, that justifies the application of the 
perturbation theory to the case of resonant coupling of the qubit to the LC circuit having 
the finite quahty factor [Qt ~ 1000). Otherwise, it is better to work out of the resonant 
conditions between the qubit and the tank. 



IV. DECOHERENCE SUPPRESSION IN A STRONGLY DRIVEN FLUX QUBIT. 

Coherent Rabi oscillations in a fiux qubit have been observed recently by Chiorescu et 
al [12]. Relaxation and dephasing times, r^eiax and r<^, of the undriven qubit have been 
measured in the process together with a decay time, TRaUi of the strongly driven system: 
Treiax = 900ns, = 20ns, and TRabi = 150ns. It is of great interest to analyze these data 
using the theoretical formulas obtained above. 
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Comparing the Hamiltonian (2) with the Hamiltonian employed in Ref. [12] wc find for 
the Rabi frequency (7): ^Ir = (A/2c<Jo)em.ui- Here luq = V + is energy sphtting of the 
qubit which is exactly equal to the circular frequency of the driving force, 27r/ : luq — 27r/, 
and Cmw is the bias-modulation amplitude linearly proportional to the amplitude of the 
microwave field [32]. 

We suppose here that decohcrcncc, dephasing and relaxation of the qubit originate from 
bias fiuctuations Q{t). These fiuctuations, which also induce fiuctuations of the qubit energy 
splitting, are characterized by the spectral density S{u!) (14) and the susceptibility 
(15). Without the external driving force the dephasing and relaxation rates, T2^eq T{'eq, 
respectively, are proportional to the spectral density S{u!) (1): 

^^^^2^:1'^^'^°^' 

T,-^,= It,-^,+ ^^S{u;c^O). (33) 

These time scales correspond to the equilibrium relaxation time Ti eg = Treiax and to the 
free-evolution dephasing time T2,eq = T<p measured by Chiorescu ct al. [12]. The pronounced 
difference between them {t~J^^ — 1.11 x 10^s~^,r^^ = 50 x 10^s~^) points to the fact that 
the slow fiuctuations of the bias give a dominant contribution (~ S{uj ^ 0)) into dephasing 
of the fiux qubit. Fortunately, the low-frequency noise does not contribute to the equilibrium 
relaxation rate 7i~eV However, it is not the case for the driven qubit. It follows from Eq.(25) 
that the nonequilibrium decay rate T^^, corresponding to the decay rate r^^^ [12], contains 
the term proportional to the spectrum of environment fluctuations at the Rabi frequency 
Sinn) : 

rr' = ^sM + j^s(si„). (34) 

Here we have taken into account that the Rabi frequency is much less than the energy 
splitting luq : Qji <^ coq (from Ref. [12] we have 0.^/271 ~ IQQM H z,ujq/ 2tx ~ A/27r ~ 
"iAGHz). The presence of low- frequency term in the decay rate (34) can explain the sub- 
stantial difference between the equilibrium relaxation rate, T^ l^ = 1.11 x 10®s~^, and the 
decay rate of the driven qubit, Tf ^ = r^^^ = 6.66 x 10^s~^. This decay rate incorporates the 
high-frequency contribution (3A^/8a;o)5'(a;o) = (3/4)r~]„3. = 0.83 x 10^s~^ as well as a term 
{e^ /AujI)S{VLii) = 5.83 x 10^s~^ which is determined by the spectrum of noise fluctuations 
taken at the Rabi frequency. 

It should be emphasized that the contribution of low-frequency fiuctuations into the 
dephasing rate of the undriven qubit 72~gg, which is proportional to the spectrum at almost 
zero frequency, S{uj ~ 0), is much higher: {e^ /2ujI)S{uj ~ 0) ~ T^^—^-^'^reiax — 49.4xl0^s~-^. 
From these data we can extract ratio between two values of the spectrum of environment 
fiuctuations: 

S{^-^) _ 2r-i - r-J^ ^ ^ ^ ^^^^ 



relax 



These estimations suggest a pronounced frequency dispersion of the spectrum S{uj) in the 
range from almost zero frequencies up to the frequencies of order ^Ir/2t: ~ IQQMHz. A 
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drastic decrease in the spectrum of environment fluctuations points to the fact that the ef- 
fective correlation time r,. of the environment, r^, is more than the period of Rabi oscillations 
(~ 10ns), Tc > 10ns. However, the correlation time Tc is probably less than 100 ns, as it fol- 
lows from the results of spin-echo experiments carried out by Chiorescu et al. [12]. It would 
be interesting to measure a dependence of the decay time TRaU on the amplitude of the exter- 
nal microwave field in the whole range of Rabi frequencies to restore a frequency dependence 
of the spectrum S{uj). For the flat spectrum, when S{Q.r) = S{uj ~ 0), the decay rate would 
be signiflcantly higher: T^ji^t = 25.3 x 10^s~^ (compared to Tj^^^ = 6.66 x 10^s~^), so that 
for the decay time of the driven qubit we would have : Tijiat = TRaUjiat = 39.5ns = TRabi/^S 
(in a comparison with the experimentally measured value Ti^bi — 150ns). The difference 
between Tjiabi and TRaujiat can be interpreted as a result of decoherence suppression (in 4 
times) by the strong microwave fleld acting on the flux qubit. The reason for this suppression 
lies in the fact that the low-frequency noise is averaged over the fast Rabi oscillations. 

Finally, we have to mention that for the strongly driven qubit it is not clear how to define 
separately the dephasing and the relaxation times. Here again we have two time scales Ti 

(25) and F"^ (19) which are related according to Eq.(29). It follows from the experimental 
data [12] that the level populations of the flux qubit are equalized for the time Ti = 150ns 

(26) , whereas establishing of steady-state polarization of the qubit in the rotating frame of 
reference (27) is determined by both time scales, Ti and F~^, with F~^ ~ 80ns. 

V. CONCLUSIONS 

We have considered a dissipative dynamics of the two-state system (the qubit) driven by 
the strong resonant fleld. The generalized Bloch equations (17) have been derived for the 
case of a weak interaction between the dynamical system and the heat bath. We have shown 
that a decay of Rabi oscillations and an establishment of the qubit polarization along the 
external oscillating field are determined by two different damping rates (25) and F^ (19) 
which are related by the simple expression (29). These results are true for weak coupling 
between the qubit and the heat bath and are not restricted to the case of Gaussian statistics 
of the heat bath fluctuations. It is demonstrated also that the Rabi oscillations eliminate 
the divergency in the decoherence rate of the biased qubit caused by the 1// noise. We 
have analyzed the necessary conditions for the observations of the Rabi oscillations in the 
phase qubit coupled to the LC-circuit having the resonant frequency which is equal to the 
Rabi frequency. It is emphasized that recent measurements of decoherence and relaxation 
times in a driven flux qubit can be regarded as an experimental evidence of decoherence 
(dephasing) suppression by the strong resonant fleld. 
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